Structural and T c inhomogeneities inherent to doping in La 2 _ x Sr x Cu04 
superconductors and their effects on the precursor diamagnetism 
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The inhomogeneities inherent to the random distribution of Sr dopants in I^-^Sr^CuCU super- 
conductors are probed by measuring the x-ray diffraction linewidths and the Meissner transition 
widths, and then consistently explained on the grounds of a simple model in which the local Sr 
content is calculated by averaging over distances close to the in-plane electronic mean free path. 
By taking into account these intrinsic bulk inhomogeneities with long characteristic lengths (much 
larger than the superconducting coherence length amplitudes), the precursor diamagnetism mea- 
sured above T c , a fingerprint of the superconducting transition own nature, is then explained for all 
doping levels in terms of the conventional Gaussian-Ginzburg-Landau approach for layered super- 
conductors. These results also suggest that the electronic inhomogeneities observed in the normal 
state by using surface probes overestimate the ones in the bulk. 

PACS numbers: 74.25.Ha, 74.40.+k, 74.62.Dh, 74.72.Dn 



I. INTRODUCTION 

It is now well established that the presence of dopants 
in high-Tc cuprate superconductors (HTSC), which so- 
deeply affect the behaviour of these materials by chang- 
ing their carrier densities, may also cause structural, elec- 
tronic, and critical temperature inhomogeneities, with 
different spatial distributions and characteristic lengths. 1 
How these intrinsic inhomogeneities inherent to doping 
affect the HTSC has now become a central open ques- 
tion of their physics^ whose interest is also enhanced 
by the existence of a wide variety of doped materials 
which may be also deeply affected by inhomogeneities 
associated with doping^ However, in spite of their rel- 
evance, aspects as important as the spatial distribution 
of these inhomogeneities or their characteristic lengths 
are still not well settled. This last is crucial, in par- 
ticular to understand how they affect the superconduct- 
ing transition, at present another open and much de- 
bated issue also related to the controversy about the na- 
ture of the pseudogap in underdoped cupratesji^ These 
lengths could be of the order of the in-plane supercon- 
ducting coherence length amplitude, £ab(0) (typically 2 to 
4 nm in La2- x Sr x Cu04), as is the case of the nanoscale 
electronic disorder observed on the surface of different 
HTSC by using scanning tunneling microscopy and spec- 
troscopy (STM/S)**™™* In this case, these intrinsic inho- 
mogeneities would deeply affect the own nature of their 
superconducting transition. In contrast, if they have in 
the bulk long characteristic lengths, much larger than 
Cab(0); they will just round the critical behaviour of any 
observable measured around T C Z 

The present debate involving the entanglement be- 
tween intrinsic inhomogeneities at different lengths 
scales, the pseudogap in the normal state and the super- 
conducting fluctuations, is particularly well illustrated 
by the diamagnetism observed above the superconduct- 
ing transition: Earlier magnetization measurements in 



an underdoped cuprate (Lai.gSro.iCuO^ already sug- 
gested that the precursor diamagnetism measured in 
the bulk is conventional, i.e., that not too close to the 
transition it could be described in terms of the conven- 
tional Gaussian-Ginzburg-Landau (GGL) approach for 
layered compounds.-^ The anomalies observed in some 
cases, in particular under low magnetic fields, could 
be easily explained in terms of T c -inhomogeneities with 
long characteristic lengths [much larger than £ a &(0)], 
which do not directly affect the superconducting tran- 
sition own nature^^ Measurements of the in-plane 
paraconductivity 11 and of the heat capacity around T c 12 
in these underdoped cuprates, together with analyses of 
the effects of T c -inhomogeneities with long characteristic 
lengths and of the background choice ^ 13 ! 14 fully agree 
with these conclusions. In particular, they confirm that 
above a reduced temperature of the order of e c = 0.5 the 
superconducting fluctuation effects vanish, as predicted 
by taking into account the limits imposed by the uncer- 
tainty principle to the shrinkage of the superconducting 
wavefunctionJ^ 

In the last years, however, different authors proposed 
that the anomalies observed in the precursor diamag- 
netism, mainly in measurements under low magnetic field 
amplitudes, provide compelling evidence that the super- 
conducting transition of underdoped cuprates is uncon- 
ventional, i.e., non describable in terms of the GGL ap- 
proach for homogeneous layered superconductorsi«^— 
Nevertheless, these authors propose very different ori- 
gins for such anomalous precursor diamagnetism, includ- 
ing the presence of superconducting phase fluctuations 
up to the pseudogap temperature or the existence in the 
bulk of intrinsic inhomogeneities with short characteristic 
lengths, as those observed by using surface probesi^^ 
These different proposals enhance the interest of a thor- 
ough study of the inhomogeneities associated with doping 
and of their effects on the precursor diamagnetism. 

In this paper, the T c -inhomogeneities inherent to the 
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random distribution of Sr dopants in I^-xSr^CuC^ su- 
perconductors are determined as a function of the doping 
level by applying to samples with high chemical quality 
two techniques that probe the bulk: powder x-ray diffrac- 
tion (XRD) and low-field magnetization^ The observed 
XRD lincwidths and T c dispersion are then explained on 
the grounds of a simple model in which the local Sr com- 
position within the CuC>2 planes results from the aver- 
age over distances of the order of the in-plane electronic 
mean free path, of the order of 30 nm. This length is 
one order of magnitude larger than both £ Q i,(0) and the 
characteristic length of the electronic disorder observed 
on the surface of different HTSC by using STM/Si^^ 
The long characteristic length of the T c -inhomogeneities 
inherent to doping is then fully confirmed by measure- 
ments of the precursor diamagnetism: The effects of these 
T c -inhomogeneities may be quenched by just applying a 
magnetic field of moderate amplitude (of the order of 10% 
the upper critical field amplitude), which shifts the aver- 
age T C (H). Then, the remaining rounding effects on the 
magnetization above T C (H) may be explained at a quan- 
titative level in terms of the conventional GGL approach 
for layered superconductors. 

II. EXPERIMENTAL 

A. Preparation of the samples and characterization 
of their structural disorder 

For our study, we have chosen the I^-xSr^CuO^ sys- 
tem because it allows to cover the entire range of doping 
by varying the Sr content. Another central experimen- 
tal aspect of our work here is to separate the chemical 
inhomogeneities inherent to doping from those extrinsic 
associated with a deficient growth process. So, to mini- 
mize these last we have decided to use powder samples 
because the Sr distribution may be very efficiently im- 
proved by successive processes of reaction and subsequent 
grinding. In addition, their small grain size (4 ± 2 /im) 
minimizes the presence of extrinsic structural inhomo- 
geneities as compared to the typical mm-sized single crys- 
tals needed to perform the magnetization experiments. 
Finally, the difficulties associated with the random ori- 
entation of the individual grains and finite-size effects 
may be easily overcome. 

The polycrystalline L^-nSr^CuO^ samples with 
0.06< x <0.275 were prepared by reacting in air at 
930° C for 20 h stoichiometric proportions of thoroughly 
mixed powders of La2C>3, SrCC>3 and CuO (99.99% pu- 
rity). The resulting polycrystalline samples were ground 
and reacted again up to ten times. A small amount 
of each composition was taken after every two grind- 
react (G-R) steps in order to monitor the sample homo- 
geneity. XRD measurements, performed with a Philips 
diffractometer equipped with a Cu anode and a Cu-K Q 
graphite monocromator, already excluded the presence 
of impurity phases after the second G-R step. An ex- 



ample of XRD pattern (corresponding to x = 0.15) is 
presented in Fig. 1(a). The resulting dependence of the 
lattice parameters on the average x value (x) is shown 
in Fig. 1(b), where the tetragonal-orthorhombic transi- 
tion at x « 0.10 is clearly seen. When analyzed in de- 
tail, the diffraction peaks present a width significantly 
larger than the instrumental one. As an example, in 
Fig. 1(c) we present a comparison between the full- width 
at half-maximum (FWHM) of some diffraction peaks for 
the x — 0.15 sample after 10 G-R steps, with the in- 
strumental FWHM. These FWHM values were obtained 
by fitting to the experimental peak profiles a pseudo- 
Voigt function [Fig. 2(a)]. It is worth mentioning that 
the linewidths decrease with the successive G-R steps: in 
the example of Fig. 2(a), after 2 G-R steps the Cu-K ai 
and K a2 components of the (008) and (400) lines over- 
lap, but after 10 steps their width is reduced and both 
components are resolved. However, as clearly shown in 
Fig. 2(b), no further reduction is observed above 8 G-R 
steps. The line broadening in excess of the instrumental 
one (A#) may be then considered as intrinsic and may be 
attributed to spatial variations of the lattice parameters 
(Act, Ab and Ac) around their average values (a, b and 
c), caused by the random distribution of Sr dopants. 

For a general orthorhombic structure, the linewidth 
A9hki (corresponding to a diffraction line with indexes 
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FIG. 1. a) Example of x-ray diffraction pattern, correspond- 
ing to the sample with x = 0.15 after 10 G-R steps, b) Result- 
ing x dependence of the lattice parameters, c) Example of the 
XRD linewidths, corresponding to the x = 0.15 sample after 
10 G-R steps, compared with the instrumental linewidths as 
determined from measurements in a reference LaB6 powder 
sample (the line is a fit of a Caglioti function). 
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h, k and I) is related to Aa, Ab and Ac through 
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where A is the x-rays wavelength. Taking into account 
the dependence of the lattice parameters on the average 
cc-value, Eq. (JTJ) may be re-written as 
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which allows to obtain the x variation about x from 
each linewidth. In Fig. 2(c), we present the x depen- 
dence of Ace, as results from the analysis of lines within 
60° < 28 < 90°, in the samples after G-R steps. We dis- 
carded lower-26* lines because of the overlapping of the 
K ai and K a2 components (see Fig. 3, left panel), and 
higher-20 lines because of their lower intensity. We also 
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FIG. 2. (color online) a) Examples (for x=0.20) of the XRD 
linewidths reduction with the successive G-R steps (differ- 
ent data sets are vertically displaced for clarity). The peaks' 
FWHM were obtained by fitting a pseudo-Voigt function 
(solid lines) taking into account the presence of both Cu-K a \ 
and K a 2 wavelengths. As shown in (b), these FWHM satu- 
rate to a value significantly larger than the instrumental one 
after 8 G-R steps, c) x dependence of the corresponding in- 
trinsic Ax. The line corresponds to Eq. ([5} evaluated with 
L = 75dcu (the shaded area represents the uncertainty in this 
value) . 
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FIG. 3. (color online) Left panel: Example of the overlapping 
of the Cu-K ai and K Q2 diffraction lines at low diffraction 
angles (note that different datasets are vertically displaced). 
Right panel: Example of how the tetragonal-orthorhombic 
transition for x < 0.10 splits diffraction lines with nonzero h 
and k indexes, making difficult the linewidths determination. 



excluded samples with x < 0.10 because the tetragonal- 
orthorhombic transition splits the lines with nonzero h 
and k indexes, introducing a large uncertainty in the A9 
determination (see Fig. 3, right panel). In the range stud- 
ied Ax ~ 0.01, in agreement with previous estimations, 27 
and it slightly increases with x. 



B. T c inhomogeneities associated with the random 
distribution of dopants, as probed by low-field 
magnetic susceptibility 

Due to the strong T c dependence on x in these 
materials^ the intrinsic compositional inhomogeneities 
measured above could lead to the existence of T c vari- 
ations (characterized by a width AT C ) around the av- 
erage T c . Here we probe these intrinsic T c inhomo- 
geneities through the temperature dependence of the 
low-field magnetic susceptibility measured under field- 
cooled (FC) conditions, X (^)j by using a commercial 
DC magnetometer (Quantum Design, model PPMS). An 
example corresponding to x = 0.08 is shown in Fig. 2] 
The magnetic field used in these measurements was al- 
ways below 1 mT which, as shown in Fig. Ufa), ensures 
that the contribution to AT C associated with the tran- 
sit through the mixed state is negligible. AT C was ob- 
tained as twice the high-temperature half-width at half- 
maximum of the dx F / dT versus T curve. In this way, 
we elude the extrinsic rounding associated with the com- 
petition between the grains size and the magnetic pen- 
etration length (which is appreciable mainly below T c ). 
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Also, as in this region \x FC \ "C 1, demagnetization ef- 
fects were neglected. As may be clearly seen in Fig.rjjb), 
AT C saturates after 8 G-R steps. This result, common to 
all compounds studied, is in agreement with the above 
XRD results, confirming the presence of an intrinsically 
inhomogeneous charge distribution in the CuO-2 planes 
associated to the random distribution of dopants. The x 
dependence of AT C /T C for the samples after 10 G-R steps 
is presented in Fig. [5l As expected from the characteris- 
tic T c {x) dependence (inset of Fig. O AT C /T C is larger 
in the most underdoped and overdoped compounds. 
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FIG. 4. a) Example (for the x = 0.08 sample after 10 G-R 
steps) of the as- measured X FC (T) under different applied 
magnetic fields. As shown in the inset, when fioH < 1 mT 
the transition width is almost unaffected by the use of a finite 
magnetic field, b) Example (for x = 0.08) of the dependence 
of the low-field X F ° (T) on the number of G-R steps. The 
procedure used to obtain T c and AT C from d\ FC /dT (solid 
lines) is also indicated. Inset: AT C dependence on the number 
of G-R steps. 



FIG. 5. x dependence of AT C /T C and of T c (inset), as ob- 
tained from the low-field \ FC '(T) measurements. The dashed 
line is a polynomial fit to data in the literature. The solid line 
corresponds to Eq. <(Sj evaluated with L — 75dc u (the shaded 
area represents the uncertainty in this value). 



III. A MODEL FOR THE INHOMOGENEITIES 

ASSOCIATED TO THE RANDOM 
DISTRIBUTION OF DOPANTS: LENGHT SCALE 
FOR THE HOLE CONCENTRATION 
INHOMOGENEITIES 

Our results on AT c {x) may be accounted for by just 
assuming that the random distribution of Sr ions in- 
duces spatial variations of the hole concentration within 
a Cu02 plane, with a characteristic length L. The lo- 
cal hole concentration p (holes/CuC^) (or, equivalently, 
the local x value) may be calculated by averaging over 
a "control" circle of radius L within the Cu02 plane, 
centered on that site, 

p = x = N St $j±. (3) 

Here N$r is the number of Sr ions within that circle, 
and dcu ~ 0.38 nm is the Cu-Cu shortest distance. 
The smaller the L value, the larger will be the statis- 
tical differences between the local and the average hole 
concentrations^ In the upper row of Fig. [5] we present 
the p (or x) distribution in a 10 3 x 10 3 d|, u area within a 
CuC>2 plane, as calculated on the basis of this model for 
several x values representative of the over-, optimal- and 
under- doped regimes. These maps were calculated by 
taking into account that there are two (La/Sr)0 layers 
associated with each CuC>2 plane, and that the proba- 
bility for a La ion to be substituted by a Sr ion is x/2. 
Also, we used L = 75<icu ~ 28 nm, which is the value 
leading to the best agreement with the above experimen- 
tal results (see below) . As expected, p varies over lengths 
of the order of L, and the spread Ap around the corre- 
sponding average value p increases with doping. As L 
is much larger than the in-plane superconducting coher- 
ence length amplitude [£ a f,(0) « 2 — 4 nm, see also below], 
and by virtue of the T c (x) relation, this p (or x) spatial 
distribution leads to the T c spatial distribution shown in 
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FIG. 6. (color online) Upper row: x (or p) distributions in a 10 3 x 10 3 dr„ area within a Cu02 plane for some representative 
x values, as calculated in terms of the model described in the text with L = 75dcu- The corresponding T c distributions are 
shown in the lower row. 



the lower row of Fig. [6] As may be clearly seen, T c also 
varies on the scale of L, but now the spread in T c val- 
ues is larger in compounds where dT c /dx is larger, i.e., 
highly underdoped or overdoped. 

We may now obtain expressions for Ap and AT C , by 
noting that the probability that iVg r out of the Nl sites 
for the Sr ions within the control circle are occupied is 
given by the binomial distribution 



N L 



an 



N L -N S 



(4) 



Here Nl may be approximated by the integer clos- 
est to 2ttL 2 /cCq , the factor 2 accounting for the two 
(La,Sr)0 layers associated with each CuC>2 layer. In 
the limit of large Nl (more precisely, if Nlx/2 3> 
5, which for L/dcu = 75 is well satisfied), Eq. (g]) 
may be approximated by a normal distribution cen- 
tered in N Sr = xirL 2 /d 2 Cn and with FWHM AN Sr w 
1.66y/N L x(l - x/2). By combining this with Eq. 
the corresponding hole density distribution is also nor- 
mal, centered in p — x and with FWHM 



Ap = Ax^ 1.33^i 



(5) 



In turn, as for the L value used in the simulations Ax <C 
x, the corresponding AT C may be approximated by 



AT c (x) w Ax(x) 



dx 



(6) 



The solid line in Fig. [5] is the best fit of Eq. © to 
the experimental data. The agreement is good, taking 
into account the experimental uncertainties, and it leads 
to L = (75 ± 25)c?cu ~ 28 ± 9 nm. It is remarkable 
that Eq. ([5]) evaluated with the same L value is in ex- 
cellent agreement with the Ax(x) dependence deduced 



from x-ray diffraction [Fig. 2(c)], which suggests a strong 
correlation between the structural and electronic disor- 
ders. Also, this L value is in reasonable agreement with 
a universal length scale recently proposed for the charge 
density inhomogeneity in these compounds for temper- 
atures around T c (~ 30 nm, see Fig. 2 in Ref. |3fj). It 
is, however, larger than the length scale for the spatial 
variations in the local density of states as determined by 
STM (5-10 nm at 4 K)J-2i 



IV^ FLUCTUATION-DIAMAGNETISM ABOVE 

T c : A CHECK OF CONSISTENCY FOR THE 
LENGTH SCALE OF THE INHOMOGENEITIES 
INHERENT TO DOPING 

One may also probe the characteristic length of 
the T c inhomogeneities associated with doping through 
measurements of the fluctuation-induced magnetization, 
Mfi, above T c . This observable is very sensitive to 
the domain size relative to the in-plane superconduct- 
ing coherence length, £ a b>^ for reduced temperatures 
e = ln(T/T c ) such that Cah(e) > L, Mft would have a 
zero-dimensional (0D) character, otherwise it would be 
two-dimensional (2D). By using the mean-field relation 
Cab(e) = £ a fc(0)£ -1 / 2 , the 2D-0D crossover would happen 
for e ps Cab(O)/-^ 2 - I n tne case that L ~ £ab(0) (as sug- 
gested by STM) all the accessible e-range would be in 
the 0D regime. 

Some examples of the M(T) curves above T c for differ- 
ent doping levels (as measured with a Quantum Design 
SQUID magnetometer) are presented in Fig. [7] From 
these measurements Mfi(T) was obtained by subtract- 
ing to M(T) the normal state contribution, which in 
turn was determined by fitting a Curie-like function 
[Mback(T) = A + BT + C/T, being A, B and C fit- 
ting constants] in a temperature interval well above T c , 
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where fluctuation effects are expected to be negligible. 
Note that whereas the background uncertainties may 
have a dramatic effect on the extracted Mfi at high re- 
duced temperatures (e > lO -1 ), their influence is small 
for e < 10" 1 , the relevant region for our present analy- 
ses. Some examples of the resulting Mfi(e) (for conve- 
nience over HT) for different doping levels are presented 
in Fig. [51 As may be clearly seen, under a 0.1 T mag- 
netic field [much smaller than Hq 2 (0), the upper criti- 
cal field for H L ab extrapolated to T = K], \Mfi\ 
presents an upturn as e —> 0. This behavior may be 
easily explained as a consequence of the T c distribution, 
and should occur below e; n h ~ 2AT c /T c f^ From the 
AT C data in Fig. El the resulting Ei n h (dashed bars in 
Fig. [5]) are in good agreement with the onset of the ob- 
served enhanced |M//|. In the presence of a finite mag- 
netic field, the region affected by inhomogeneities, i.e., 
around T C (H), is shifted to lower temperatures. This 
allows to compare the theoretical results for homoge- 
neous materials with the experimental data in a wider 
temperature region above T c (0). Quantitatively, it may 
be approximated e inh (H) ~ 2AT C /T C - H/H£ 2 (0)H As 
H H^ 2 (0) ~ 20 - 80 T in these materials^ a 5 T mag- 
netic field will allow to study Mfi free from the effect of 
T c inhomogeneities down to e ~ for 0.1 < x < 0.2. 

In presence of a finite magnetic field perpendicular to 
the CuC>2 planes, Mjj in the framework of the 2D Gaus- 
sian Ginzburg-Landau (GGL) approach is given by^ 



M fi = 
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Here T and ip are the gamma and digamma functions, 
h = H/Hg 2 (0) the reduced magnetic field, s « 6.6 A the 
periodicity length of the superconducting layers, / the 



60 



x=0.20 




100 



T(K) 



FIG. 7. Some examples of the as-measured M(T) in the 
normal state. The effect of superconducting fluctuations is 
clearly seen as a rounding close to the corresponding T c (indi- 
cated by arrows). The normal-state background (solid lines) 
was determined by fitting a Curie-like function above 60 K. 
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FIG. 8. Reduced-temperature dependence of Mfi (over HT) 
for different doping levels. The dramatic effects of T c inho- 
mogeneities inherent to doping manifested at low fields are 
quenched by the application of a 5 T magnetic field, allowing 
a direct comparison with the 2D-GGL theory (lines). Dashed 
bars are the expected onset of the \Mfi\ upturn due to these 
T c inhomogeneities at low fields. Inset: The resulting £ab(0) 
values are presented against T c _1 . 



effective superconducting fraction,-- 4 ks the Boltzmann 
constant, (f>o the flux quantum, and e c « 0.5 the total- 
energy cutoff constant i 1 ^ In the low magnetic-field limit 
(h -C e) Eq. |J7J| may be approximated by 



60o s 



1 



(8) 



which in absence of a cut-off (e c — » oo) reduces to the 
well known Schmid result for 2D materials.— As in these 
materials the anisotropy factor, 7, is in the range 10-25, ?6 
the fluctuation magnetization when H \\ ab (of the order 
of ~ Mji/j 2 , see e.g. Ref.E3) is negligible. Then, M fl in 
a powder sample with the grains randomly oriented may 
be approximated by the angular average^ 



(Mfi) = \ 



M f \{Hcos6)cos9sm6d0, (9) 



where 9 is the angle between the c-axis and H . As may be 
clearly seen in Fig.[8j the comparison of the 2D-GGL the- 
ory (solid lines) with the 5 T data is excellent down to the 
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lowest accessible reduced temperature (e ~ 1CP 2 ), the re- 
sulting £ afc (0) = [(/) /27r/x iJp 2 (0)] 1/2 being about 2-4 nm 
(inset in Fig. HJ in agreement with recent estimates.— 
The 2D nature of Mfi is consistent with a length scale 
for the in-plane modulation of the hole-concentration well 
above the nanometer scale, suggesting that the STM 
resultsii 2 -'^— could not be representative of the bulk. 

V. CONCLUSIONS 

The results we have obtained in the prototypical 
La2- 2: Sr 2 ;Cu04 system provide unambiguous answers to 
some of the questions we have stressed in the intro- 
duction: i) The unavoidable inhomogeneities when dop- 
ing the cuprate superconductors have long characteris- 
tic lengths [much larger than £ a fe(0)] and they are uni- 
formly distributed in the bulk, ii) These intrinsic in- 
homogeneities may deeply affect the bulk magnetization 



measurements around T c , but once they are taken into 
account the diamagnetic transition may be described for 
all doping regions in terms of the Ginzburg-Landau ap- 
proach for layered superconductors, i.e., from a phe- 
nomenological point of view the intrinsic diamagnetic 
transition in the HTSC is conventional. Additionally, 
our results do not support the presence of fluctuating vor- 
tices well above the superconducting transition, even in 
the pseudogap regime, and they also suggest that the dif- 
ferent inhomogeneities observed by using surface probes 
overestimates the ones in the bulk. 
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